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摘 要 : 基于 2000 一 2019 年 的 林 草 植被 覆盖 度数 据 分 析 了 河套 灌区 近 20 a 的 林 草 植被 时 空 变化 特征 ， 


利用 地 


理 探测 需 定 量 分 析 了 各 驱动 因子 对 林 草 植被 覆盖 的 影响 。 结 果 表 明 :(1) 近 20a 河 套 灌区 林 草 地 面积 总 体 呈 减少 
态势 ,减少 量 为 966.15 km , 林 草 植被 覆盖 度 等 级 整体 呈 上 升 趋势 ,多 年 平均 林 草 植被 覆盖 度 为 28.3% ,空间 尺 度 上 ， 
河套 灌区 林 草 植被 覆盖 度 呈 现 出 东北 向 西南 递减 的 空间 变化 特征 ;(2) 近 20 a 低 林 草 植被 覆盖 度 .中 林 草 植被 覆盖 


度 等 重心 主要 位 于 永 济 灌 域 东北 侧 和 义 长 灌 域 西南 侧 , 各 级 林 草 植被 覆盖 重心 表现 为 从 东北 向 西南 方向 递减 ， 


均 呈 现 出 向 西南 方向 迁移 的 态势 ,各 时 段 迁 移 距 离 介 于 0.71~15.46 km; (3) 距 水 域 距离 地 下 水 埋 深 .气温 .降水 是 
影响 区 域内 林 草 植被 覆盖 度 的 主导 环境 因子 , 且 解 释 力 分 别 为 0.427、0.439、0.318、0.368; 距 水 域 距离 .地 下 水 埋 深 、 


年 平均 降水 与 其 他 因子 的 交互 作用 普遍 较 高 ,区 域 含水 量 的 大 小 是 影响 河套 灌区 林 草 地 生长 状况 的 主要 驱动 力 ， 
研究 结果 可 为 灌区 自然 植被 的 保护 与 恢复 提供 相应 的 科学 依据 。 


关键 词 : 河套 灌区 ; 林 草 植被 覆盖 度 ; 时 空 变化 ; 驱动 


林 草 作为 陆地 生态 系统 的 重要 组 成 部 分 之 一 ， 
在 全 球 物质 循环 及 能 量 流动 中 发 挥 着 重要 的 作用 ， 
对 缓解 地 区 水 土 流失 具有 明显 的 生态 效益 "”"。 植 
被 覆盖 度 (Fractional Vegetation Cover, FVC) 可 以 作 
为 衡量 植被 生长 状况 的 一 个 重要 的 指标 ,能 够 反映 
出 区 域 自然 生态 系统 中 存在 的 植被 退化 等 各 种 环 
境 问题 ,是 估算 和 监测 区 域 生 态 系统 稳定 性 和 功能 
的 重要 参数 ,很 多 学 者 基于 植被 覆盖 度数 据 应 用 
不 同方 法 分 析 了 区 域 植被 的 生长 变化 情况 及 驱动 
力 分 析 , 大 量 人 研究 发 现 气 象 因 素 ( 气 温 、 降 水 ) 对 植 
被 生长 变化 的 解释 力 较 强 ”中 。 如 马 晓 妮 等 "通过 
对 黄土 高 原 的 而 砂 岩 区 植被 覆盖 度 环境 驱动 因子 
的 量化 分 析 表 明 ,降水 .气温 .土壤 水 分 是 影响 看 砂 
岩 区 植被 覆盖 空间 分 布 的 主导 环境 因子 。 近 年 来 ， 
更 多 的 学 者 定量 分 析 和 评价 了 自然 环境 因子 对 植 
被 生长 的 影响 ,而 对 植被 生长 与 人 类 活动 的 响应 关 
系 探究 较 少 ,实际 上 ,人 类 活动 对 植被 生长 与 变化 
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因子 ; 地 理 探 测 器 


DEADEN 

地 理 探 测 器 是 基于 一 组 统计 学 方法 来 探测 空间 
分 异性 ,揭示 其 背后 驱动 力 的 空间 统计 方法 ““ 呈 。 
祝 聪 等 光 通 过 地 理 探 测 需 对 2006 一 2016 年 岷江 上 
游 植被 覆盖 度 的 时 空 变化 进行 了 驱动 力 分 析 , 发 现 
岷江 上 游 植被 覆盖 度 主 要 受 海 拔 .气温 .土壤 类 型 、 
降水 4 个 因子 的 影响 。 许 多 研究 对 区 域 植被 覆盖 度 
时 空 变化 的 驱动 力 分 析 有 着 重要 的 意义 ,但 大 多 只 
是 对 驱动 因子 对 植被 覆盖 度 的 解释 力 进 行 了 定量 
分 析 , 缺 乏 对 影响 植被 生长 的 驱动 因子 适宜 条 件 的 

河套 灌区 大 部 分 区 域 处 于 蕊 漠 化 边缘 地 带 , 导 致 
该 地 区 的 草地 、 林 地 等 自然 植被 处 于 退化 的 边缘 ™， 
并 且 作 为 我 国 较 大 的 农 牧 交 错 带 ,常年 的 耕作 和 灌 
溉 ,导致 其 地 下 水 位 下 降 和 土壤 盐 碱 化 加 剧 ,从 而 
加 重 了 生态 系统 的 复杂 性 和 脆弱 性 , 若 不 子 以 有 效 
保护 , 极 易 向 荒漠 化 演变 2 。 近 年 来 ,不 少 学 者 对 
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包含 河套 灌区 在 内 的 黄河 流域 和 蒙 陕 甘地 区 的 植 
被 覆盖 进行 了 研究 ,其 中 ,张志强 等 ”分 析 了 


东西 跨度 约 为 250~260 km, 南北 跨 度 约 为 45~60 
km。 江 区 地 形 平坦 ,整体 海拔 介 于 911~1167 m 之 


2000 一 2019 年 黄河 流域 植被 覆盖 度 时 空 分 布 格局 ， 
对 制定 黄河 流域 生态 保护 方案 提供 了 数据 支撑 ; 苏 
丽 德 等 的 对 1986 一 2019 年 乌拉 特 后 旗 荒 漠 化 动态 
变化 进行 了 分 析 ,发现 了 乌拉 特 后 旗 90% 的 区 域 都 
处 于 荒漠 化 状态 。 然 而 ,以 往 的 植被 覆盖 度 的 研究 
一 般 都 是 对 区 域 总 体 植被 进行 了 分 析 ,而 对 于 以 农 
业 为 主 的 河套 灌区 ,区 域 总 体 植 被 覆盖 度 的 结论 很 
难 反映 出 自然 植被 的 生长 情况 。 因 此 本 文 以 此 为 
切入 点 ,开展 河套 灌区 林 草 植被 年 际 动 态 变 化 的 
研究 。 


1 材料 与 方法 


1.1 研究 区 概况 

河套 灌区 位 于 我 国内 蒙古 自治 区 巴 彦 深 尔 市 ， 
是 我 国 三 大 灌区 之 一 ,也 是 亚洲 最 大 的 灌区 (图 
1)。 灌 区 属 温带 大 陆 性 季风 气候 ,年 降水 量 较 少 ， 
年 平均 降水 量 约 为 135~225 mm, 主要 集中 于 7 一 9 
月 。 蒸 发 量 比较 大 ,年 平均 蒸发 量 约 为 1750 mm, 
多 年 平均 气温 约 为 5.6 'C ,灌区 北 依 阴 山 山 脉 , 南 临 
黄河 , 东 接 内 蒙古 包头 市 郊 , 西 靠 乌 兰 布 和 沙漠 。 
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间 。 灌 区 年 引 黄 水 量 约 S0x10"ms , 占 黄河 过 境 水 量 
的 六 分 之 一 ,是 灌区 重要 的 水 资源 来 源 。 
1.2 数据 来 源 及 预 处 理 

本 文 共 使 用 与 研究 相关 的 12 种 基础 数据 。 包 
括 用 于 提取 植被 覆盖 度 的 遥感 数据 ,提取 林 草 植被 
的 土地 利用 矢量 数据 以 及 用 于 分 析 的 地 形 数据 、 气 
象 数 据 、 地 下 水 数据 .CDP 和 人口 数据 、 距 水 域 . 公 
路 .居民 点 .耕地 距离 数据 ( 表 1)。 
1.3 研究 方法 
1.3.1 林 草 植被 覆盖 度 “ 利 用 像 元 二 分 模型 估算 植 
被 覆盖 度 汪 ,计算 公式 如 下 : 

NDVI-NDVI,,, 

NDVI,,, -NDVI 
式 中 :FVC 表 示 植 被 覆盖 度 ,本 文 参考 李 苗 苗 等 59 
提出 的 NDVL 和 NDVL 的 取 值 方法 ,在 NDVI 的 频 
率 累 计 表 上 取 给 定 置信 区 间 内 的 最 大 值 和 最 小 值 ， 
分 别 作为 NDVL。. 和 NDVILs, 因 此 本 文 结合 河套 灌区 
的 实际 情况 与 前 人 研究 采用 5%~95% 的 置信 区 间 佑 
算 植 被 覆盖 度 ,通过 林 草 地 矢量 数据 提取 河套 灌区 
林 草 地 植被 覆盖 度 。 

参考 中 华人 民 共 和 国 水 利 部 2008 年 颁发 的 《 土 


FVC= 


(1) 
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图 1 研究 区 概况 示意 图 
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Fig. 1 Survey map of the study area 
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表 1 双 因 子 交互 作用 结果 分 类 


Tab.1 Two-factor interaction results classification table 


数据 类 型 数据 来 源 产品 分 辩 率 (时 间 范 围 ) 预 处 理 
遥感 数据 Google Earth Engine 遥感 云 平 台 Landsat Surface Reflec- 30 m(2000—20194F ) 据 集 都 经 过 几何 校正 ,大气 校正 等 


(https://earthengine.google.com) 


土地 利用 类 内 蒙古 土地 调查 规划 院 


tance Tier 1 系列 数据 


巴 彦 淖尔 市 各 个 旗 县 的 2000 年 .2009 年 .2019 年 
型 矢量 数据 土地 利用 类 型 矢量 数据 


选取 每 年 7 一 10 月 植被 生长 旺盛 季 作 
数据 源 

上 用 ArcGIS 软件 进行 拼接 裁剪 等 处 

理 得 到 河套 灌区 土地 利用 矢量 数据 


m 
预 处 理 , 同 时 结合 区 域 的 物候 规律 ， 
选 
为 


地 形 数据 ”地 理 空间 数据 云 平 台 (http://www，ASTER GDEM 数字 高 程 30m (2009 年 .2019 年 ) iit ENVIAR AR EFT CHE Bin 、 投 影 、 
gscloud.cn ) 产品 裁剪 等 处 理 得 到 河套 灌区 高 程 数 据 

气象 数据 ” 中国 气 象 科学 数据 共享 网 站 (http: 降水 量 .气温 2009 年 .2019 年 基于 反 距 离 权重 内 插 法 生成 150 m 分 
//data.cma.gov.cn/ ) 辩 率 的 栅 格 数据 

地 下 水 数据 ”通过 自动 水 位 监测 系统 采集 ,地 下 一 2009 年 .2019 年 通过 数据 整理 与 插值 计算 得 到 河套 
水 位 观测 井 共 176 眼 灌区 150 m 分 状 率 的 地 下 水 位 栅 格 图 

GDP 和 人 口 巴 彦 淖尔 市 统计 局 巴 彦 淖尔 市 社会 统计 年 2009 年 .2019 年 基于 反 距离 权重 内 插 法 生成 150 m 分 

数据 鉴 辩 率 的 栅 格 数据 

HE Kk Sak . ZS - 一 2009 年 .2019 年 提取 土地 利用 矢量 数据 中 各 要 素 , 基 

路 、 居 民 点 、 于 欧 氏 距离 计算 每 个 像 元 到 各 要 素 

耕地 距离 之 间 的 距离 ,通过 重 采样 生成 150 m 

分 辨 率 的 栅 格 数据 


培 侵 蚀 分 类 分 级 标准 和 "与 植被 覆盖 度 分 级 要 求 的 
相关 文献 结合 野外 实际 调查 情况 ,将 植被 覆盖 度 通 
过 密度 分 割 划分 成 5 个 等 级 ( 表 2)。 

表 2 植被 覆盖 度 分 级 


Tab. 2 Vegetation coverage classification 


植被 覆盖 度 等 级 FVC 
低 植被 覆盖 度 <15% 
较 低 植被 覆盖 度 15%~30% 
中 植被 覆盖 度 30% ~45% 
较 高 植被 覆盖 度 45%~60% 
高 植被 覆盖 度 >60% 


1.3.2 FVC 变化 趋势 分 析 法 “年 均 FVC 变 化 趋势 表 
达 式 如 下 : 


n 


n n 
"Si, Sis, 
i=l 


iit (2) 
Dy by 
SUF: slope 表示 林 草 植被 覆盖 度 变 化 趋势 斜率 ;n 表 
示人 研究 时 间 序 列 的 长 度 ;变量 i 表示 年 份 序号 ;x; 表 
示 第 i 年 的 植被 覆盖 度 。 其 中 , 当 slope 值 为 正 时 , 表 
明 在 研究 时 段 内 的 植被 覆盖 度 呈 增加 趋势 ; 耕 slope 
值 为 负 值 , 则 表明 植被 覆盖 度 在 研究 时 段 内 呈 减 少 
趋势 。 
1.3.3 重心 迁移 模型 ”重心 迁移 模型 是 较为 常见 的 
分 析 模 型 之 一 , 它 能 够 直观 地 表现 出 某 个 要 素 在 空 


slope = 


间 上 迁移 的 规律 以 及 集聚 的 情况 所 。 重 心 迁 移 模 
型 计算 公式 如 下 : 


S(F.X 4) 
X,= =4___ (3) 
yes 
(FXy,) 
y= (4) 


>F. 

SUP X YRI E EAE mi E Eob ZS A FE As 
PR; FRIG i ER at 2 ROUTE s FAE xiy: 
为 第 * 年 在 第 ;个 格 网 上 植被 覆盖 度 的 重心 坐标 。 
1.3.4 FVC 驱动 因子 定量 分 析 

(1) 因子 探测 

探测 植被 覆盖 度 的 空间 分 异性 和 各 个 驱动 因 
子 对 林 草 植被 柳 盖 度 空 间 分 异性 的 解释 程度 ,用 gq 
值 度量 ,计算 公式 如 下 : 


SSW= SN,o?, SST= No? (6) 


h=1 


式 中 :为 林 草 植被 覆盖 度 或 驱动 因子 的 分 类 或 分 
区 ;WV 和 NN 分别 为 第 h 分 区 与 全 区 的 单元 数 ; o 和 
中 分 别 为 第 h 分 区 与 全 区 y 值 的 方差 。SSW 和 SST 
分 别 为 屋内 方差 之 和 与 全 区 总 方差 。g 的 值 域 为 
[0,1],g 值 越 大 则 表示 驱动 因子 x 对 植被 覆盖 度 y 


(5) 
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的 解释 力 就 越 强 , 反 之 则 越 弱 。 

(2) 交互 探测 需 

交互 探测 需 用 于 识别 不 同 驱 动因 子 之 间 的 交 
互 作 用 , 即 探测 驱动 因子 入 和 驱动 因子 马 共 同 作用 
时 对 植被 覆盖 度 的 解释 力 是 增强 还 是 减弱 。 不 同 
因子 交互 作用 的 差异 主要 体现 出 5 种 结果 ,4 
X2)<Min | q(%), (X2) | 交互 作用 非 线性 增强 ,Min[y 
(X), (X) }<q(XiNX2)<Max | q(X1), (Xo) 128 AEA A 
FAI SEZ PES aR, g(XiX2)>Max Lq(X1), qX) 26 ENE 
FRA SE TIIR , (XNXX HA EA Ah 
SEC EAE FH , (XNXX HX VE AF ERE 
增强 。 

(3) 风险 区 探测 

风险 区 探测 用 于 判断 各 驱动 因子 的 子 区 域 间 
的 属性 均值 是 否 有 显著 的 差别 ,并 且 可 以 找 出 林 草 
植被 适宜 生长 的 范围 。 风 险 探 测 用 :统计 量 来 检 
验 ,表达 公式 如 下 : 

i= Ver Ves I (7) 
=i -1) p Var(Y, a| 


N,- Ni =2 


式 中 :了 为 子 区 域 凡 内 林 草 植被 覆盖 度 的 线性 回归 
系数 ;nm 为 子 区 域 h 内 样本 的 数量 。 


2 结果 与 分 析 


21 河套 灌区 林 草 FVC 空间 分 布 特征 

图 2 为 2000 一 2019 年 5 个 时 期 的 植被 覆盖 度 等 
级 空间 分 布 图 。 河 套 灌区 在 每 年 林 草 长 势 最 旺盛 
的 7 一 9 月 林 草 植被 覆盖 度 总 体 呈 东 高 西 低 的 态势 
分 布 , 多 年 来 43.2%~51.6% 的 低 植 被 覆盖 度 林 草 主 
要 分 布 于 西部 的 乌 兰 布 和 灌 域 ,53.1%~57.4% 的 较 
高 和 高 植被 覆盖 度 林 草 主要 分 布 于 东部 的 乌 梁 素 
海 的 周边 地 区 。 这 是 由 于 西部 林 草 多 以 沙漠 一 草 
甸 的 形式 存在 ,地 下 水 位 较 低 、 稳 定性 较 弱 ,导致 林 
草 的 长 势 较 差 ,而 东部 由 于 存在 着 乌 梁 素 海 等 较为 
丰富 的 水 源 ,并 且 近 年 来 受 " 一 湖 两 海 " 政 策 的 影 
响 ,生态 保护 力度 加 强 。 但 20 a 来 西南 地 区 的 较 高 
植被 覆盖 度 和 高 植被 覆盖 度 的 数量 增加 , 而 东北 地 
区 却 有 所 减少 。 

为 了 进一步 探究 和 分 析 河 套 灌 区 林 草 履 被 变 
化 的 情况 ,对 2000 一 2019 年 不 同等 级 的 林 草 植被 覆 
盖 度 进行 等 级 转移 和 矩阵 运算 。 表 3 为 2000 一 2019 
年 河套 灌区 林 草 植被 覆盖 度 等 级 转移 矩阵 。 可 以 


看 出 ,研究 时 段 内 共有 1647.79 km 的 其 他 土地 利用 
类 型 转变 为 林 草 ,其 中 转 入 较 低 植被 覆盖 度 和 中 植 
被 覆盖 度 的 较 多 ,分别 为 652.75 km 、413 kms FEA 
2613.94 km 的 林 草 地 转变 其 他 土地 利用 类 型 ,其 中 
低 植被 覆盖 度 和 较 低 植被 覆盖 度 转 出 的 较 多 ,分别 
为 984.38 km2 1078.26 km ,并 且 转 出 部 分 主要 分 布 
于 乌 兰 布 和 灌 域 的 荒漠 地 带 。20 a 间 , 林 草地 数量 
整体 呈 减 少 的 态势 ,但 从 植被 覆盖 度 等 级 转化 来 
看 , 低 等 级 植被 覆盖 度 转 为 高 等 级 植被 覆盖 度 的 数 
量 为 520.62 km’, 高 等 级 植被 覆盖 度 转 为 低 等 级 植 
被 覆盖 度 的 数量 为 151.51 km ,由 此 看 出 ,河套 灌区 
林 草 地 生长 状况 整体 上 得 到 改善 。 
2.2 河套 灌区 FVC 时 序 变 化 特征 

如 图 3 所 示 ,河套 灌区 林 草 植被 覆盖 以 较 低 植 
被 覆盖 为 主 , 占 林 草 植 被 总 面积 的 37.01%~43.53%; 
其 次 是 低 植 被 覆盖 度 、 中 植被 覆盖 和 较 高 植被 覆 
盖 ,而 高 植被 覆盖 的 林 草 地 则 分 布 最 少 , 占 林 草 植 
被 总 面积 的 3.8%~6.24%。 其 中 在 2000 一 2009 年 低 
植被 覆盖 和 较 低 植被 覆盖 占 比 总 体 呈 下 降 态 势 , 而 
中 植被 覆盖 、 较 高 植被 覆盖 和 高 植被 覆盖 呈 上 升 态 
热 , 这 说 明 该 时 期 河套 灌区 林 草 地 的 生长 状况 得 到 
了 了 改善。 在 2009 一 2014 年 低 植 被 覆盖 和 较 低 植被 
履 盖 占 比 总 体 呈 上 升 态 势 ,而 中 植被 覆盖 、 较 高 植 
被 覆盖 和 高 植被 覆盖 呈 下 降 的 态势 ,说 明 这 时 期 林 
草 生长 状况 有 所 恶化 。 在 2014 一 2019 年 低 植被 覆 
盖 和 较 低 植被 覆盖 占 比 总 体 又 呈 下 降 态 势 , 而 中 植 
被 覆盖 、 较 高 植被 覆盖 和 高 植被 覆盖 呈 上 升 的 态 
势 ,说 明 此 时 段 河套 灌区 林 草 生长 状况 得 到 改善 。 

由 图 4 可 以 看 出 ,河套 灌区 多 年 平均 FVC 为 
0.283 ,平均 FVC 在 2002 年 .2007 年 .2009 年 .2019 
年 出 现 峰值 ,分别 为 0.337 0.323 .0.353 .0.304, 谷 值 
出 现在 2001 年 .2006 年 .2008 年 .2016 年 ,分别 为 
0.234 年 .0.278 年 .0.250 年 .0.265 年。 在 整个 研究 
时 段 内 河套 灌区 平均 FVC 呈现 出 较为 明显 的 波动 
震荡 ,平均 FVC 的 高 低 交 替 比 较 频 繁 ,多 年 平均 
FVC 呈现 为 较 低 植被 覆盖 度 ,FVC 的 线性 趋势 不 显 
著 , 整 体 上 变化 不 明显 。 相 对 变化 率 在 2001 一 2002 
年 .2005 一 2006 年 .2008 一 2009 年 .2016 一 2017 年 增 
长 的 速率 较 快 ,说 明 这 4 个 时 期 林 草 覆盖 增多 ,得 到 
了 较 好 的 改善 。 而 在 2002—2003 年 .2004 一 2005 
Œ 2009—2010 4 2014—2015 年 减少 的 速率 较 快 ， 
说 明 这 4 个 时 期 林 草 植被 覆盖 度 相 对 上 一 年 有 所 
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图 2 林 草 植被 覆盖 度 分 布 


Fig. 2 Vegetation coverage distribution of forest and grass 


#3 2000—2019 年 河套 灌区 林 草 植被 覆盖 度 等 级 转移 矩阵 


Tab.3 Grade transfer matrix of forest-grass vegetation coverage in Hetao Irrigation District from 2000 to 2019 /km 


FVC 等 级 低 较 低 中 较 高 高 转 入 
低 111.91 58.24 7.06 1.5 0.66 342.49 
较 低 182.67 237.46 44.39 8.25 3.05 652.75 
中 71.14 153.32 55.92 14.68 6.42 413 
较 高 15.77 37.2 23.01 10.54 7.26 143.93 
高 6.6 12.7 10.28 7.93 11.44 95.62 
转 出 984.38 1007.26 401.12 133.27 87.91 
减少 。 5 可 以 看 出 , 低 覆 盖 类 林 草 重心 在 整个 研究 时 段 均 
2.3 各 级 林 草 FVC 重心 迁 移 位 于 永 济 灌 域 内 ,在 2000 一 2003 年 , 低 植 被 覆盖 类 


重心 迁移 模型 主要 体现 在 方向 和 距离 特征 上 林 草 的 重心 迁移 的 距离 最 大 ,为 12.154 km; BOA 
的 变化 规律 ,可 以 反映 出 景观 时 空 变化 情况 。 由 图 。 盖 类 林 草 重心 在 整个 研究 时 有 段 均 位 于 永 济 灌 域 西 
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图 4 河套 灌区 林 草 植被 平均 FVC 变化 特征 


Fig. 4 Average FVC variation characteristics of forest and grass vegetation in Hetao Irrigation District 


南 侧 ,总 体位 于 低 覆 盖 类 的 东北 侧 ,在 2013 一 2016 
年 , 较 低 植被 覆盖 类 林 草 的 重心 迁移 的 距离 最 大 ， 
为 11.857 km; 中 覆盖 类 林 草 重心 分 布 于 永 济 灌 域 和 
义 长 灌 域 的 交界 地 带 , 总 体位 于 低 覆 盖 类 和 较 低 覆 
盖 类 的 东北 侧 ,在 2003 一 2006 年 ,中 植被 覆盖 类 林 
草 的 重心 迁移 的 趾 离 最 大 ,为 8.137 km; 较 高 覆盖 类 
林 草 重心 分 布 于 永 济 灌 域 和 义 长 灌 域 的 交界 地 带 ， 
在 2003 一 2006 年 , 较 高 植被 覆盖 类 林 草 的 重心 迁移 
的 距离 最 大 ,为 15.455 km; 高 覆盖 类 林 草 重心 分 布 于 


义 长 灌 域 的 西南 侧 , 在 2003 一 2006 年 ,高 植被 覆盖 类 

林 草 的 重心 迁移 的 距离 最 大 ,为 9.707 km。 各 级 林 草 

植被 覆盖 重心 总 体 呈 向 西南 迁移 的 态势 ,并 各 级 林 

草 植 被 覆盖 重心 表现 为 从 东北 向 西南 方向 递减 。 

24 林 草 植被 覆盖 度 变 化 驱动 力量 化 

2.4.1 探测 因子 单 因素 分 析 “为 探究 不 同 环境 因子 

与 人 类 活动 对 河套 灌区 林 草 植被 覆盖 度 空间 分 异 
4 影响 ,选取 了 9 种 探测 因子 进行 分 析 ,得 到 单 因 子 

影响 的 解释 力 4 值 (图 6) , 林 草 植被 覆盖 度 空 间 分 异 
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图 5 各 级 植被 覆盖 度 重 心 迁 移 
Fig. 5 Gravity center migration map of vegetation coverage at all levels 
各 因子 的 解释 力 排序 依次 为 :地 下 水 (0.439)> 距 水 242 探测 因子 交互 分 析 通过 交互 探测 器 分 析 了 


La iB BS (0.427 ) > KEK (0.368 ) > 气温 (0.318)> 距 耕地 


河套 灌区 的 林 草 植被 覆盖 度 的 各 个 影响 因子 间 的 


距离 (0.241)> 距 道路 距离 (0.220)> 距 居民 点 距离 
(0.174)>GDP(0.107)> 高 程 (0.026)。 其 中 ,地 下 水 、 
距 水 域 距离 .降水 .气温 的 g 值 均 大 于 0.3, 是 河套 灌 
区 林 草 植被 覆盖 空间 分 布 的 主导 环境 因子 ,而 距 耕 
地 距离 . 距 道 路 距离 .中 居民 点 距离 .GDP 的 gq 值 小 
于 0.3, 且 对 灌区 林 草 植被 覆盖 空间 分 布 的 解释 力 相 
对 较 小 ,由 此 说 明 影 响 河套 灌区 的 林 草 植被 覆盖 空 
间 分 布 主要 为 自然 因素 。 


交互 作用 ,结果 如 图 7 所 示 , 在 双 因 子 的 交互 作用 下 
会 加 强 对 林 草 植被 覆盖 度 的 解释 力 ,以 双 因 子 增强 
为 主 。 河 套 灌区 双 因 子 之 间 交 互 作用 对 林 草 植被 
覆盖 度 的 影响 较 大 的 有 距 耕 地 距离 门 地 下 水 埋 深 
(0.693)> 距 水 体 距 离 门 距 居 民 点 距离 (0.610)> 距 居 
民 点 距离 站 年 平均 气温 (0.524)> 距 道路 距离 门 距 水 
体 距 离 (0.519)= 距 水 体 距 离 门 地 下 水 埋 深 (0.519)= 
距 居 民 点 距离 门 地 下 水 埋 深 (0.519) ,其 余 因 子 交互 
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图 6 河套 灌区 林 草 FVC 空 间 分 异 驱动 因子 解释 力 


Fig. 6 Explanatory power statistics of driving factors of forest-grass FVC spatial differentiation in Hetao Irrigation District 
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图 7 2000 一 2019 年 河套 灌区 各 驱动 因子 交互 探测 解释 力 (g) 
Fig.7 Interactivity detection explanatory power (q) of driving factors in Hetao Irrigation District from 2000 to 2019 


作用 对 林 草 植被 覆盖 度 的 解释 力 均 小 于 0.5, 距 水 域 
距离 .地 下 水 埋 深 \ 年 平均 降水 与 其 他 因子 的 交互 
作用 普遍 较 高 ,表明 了 区 域 仿 水 量 的 大 小 是 影响 河 
套 灌 区 的 林 草 地 生长 情况 好 坏 的 主要 因子 。 

2.4.3 风险 区 探测 ”由 于 本 文 所 选取 驱动 因子 中 的 
GDP 和 高 程 的 单 因素 分 析 与 交互 分 析 的 解释 力 均 
比较 小 ,所 以 风险 探测 仅 列 出 了 其 他 7 个 主要 因 


素 。 如 表 4 所 示 ,河套 灌区 林 草 植被 在 147~157 mm 
的 降水 区 间 内 植被 覆盖 度 均值 最 高 ,并且 总 体 上 林 
草 植被 覆盖 随 着 降水 的 增多 而 逐渐 增 大 ,与 降水 量 
呈正 相关 关系 ,这 也 表明 河套 灌区 降水 量 的 多 少 是 
林 草 地 的 生长 状况 好 坏 的 关键 因素 。 林 草 植被 在 
8.9~9.2 C 年 平均 气温 区 间 内 的 植被 覆盖 度 均值 最 
高 , 林 草 地 植被 在 1.8~2.8 m 地 下 水 区 间 内 的 植被 覆 
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Tab. 4 Optimum range of driving factors 


探测 因子 指标 最 大 FVC 均 值 最 适 范 围 
XI 距 道 路 距离 0.387 10014~13121 m 
X2 距 耕 地 距离 0.467 1566~2254 m 
X3 距 水 体 距 离 0.327 0~670 m 
X4 距 居 民 点 距离 0.402 7537~11250 m 
X5 地 下 水 埋 深 0.499 1.8~2.8 m 
X6 年 平均 气温 0.351 8.9~9.2 °C 
XT 年 平均 降水 0.665 147~157 mm 


盖 度 均值 最 高 , 林 草 植被 在 7537~11250 m HEJS R A 
距离 区 间 及 10014~13121 m 距 道路 距离 区 间 内 林 草 
植被 覆盖 度 均 值 最 高 ,并 且 总 体 上 林 草 地 植被 的 生 
长 状况 与 距 居 民 点 距离 与 距 道 路 距离 呈正 相关 关 
系 ,表明 林 草 植被 的 生长 状况 与 人 类 活动 呈 负 相关 
关系 ,距离 人 类 活动 较为 活跃 的 地 方 越 远 , 林 草 的 
生长 状况 越 好 , 林 草 植被 在 1566~2254 m EPAF HEE 
离 区 间 林 草 植被 覆盖 度 均 值 最 高 , 林 草 植被 在 0~ 
670 m 路 水域 距离 区 间 林 草 植 被 覆盖 的 均值 最 高 。 


3 讨论 

通过 对 河套 灌区 林 草 地 长 时 序 的 探究 发 现 : 在 
研究 时 段 内 河套 灌区 的 低 植被 覆盖 度 林 草 主要 分 
布 于 乌 兰 布 和 灌 域 ,而 较 高 和 高 植被 覆盖 度 林 草 主 


年 均 温 变化 趋势 
= -0.04~0.02 mm 0.06~0.08 m 0.14~0.18 
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要 分 布 于 东部 的 乌 梁 素 海 的 周边 地 区 ,这 是 因为 西 
部 林 草 多 以 沙漠 一 草 甸 的 形式 存在 ,地 下 水 位 较 
低 .稳定 性 较 弱 ,导致 林 草 的 长 势 较 差 ,而 东部 由 于 
存在 着 乌 梁 素 海 等 较为 丰富 的 水 源 , 并 且 受 "一 湖 
两 海 ” 等 政策 的 影响 ,生态 保护 力度 加 强 , 林 草 植被 
窗 盖 度 总 体 呈 东 高 西 低 的 态势 分 布 ,并 且 多 年 河套 
灌区 的 林 草 植被 覆盖 度 等 级 呈 s 上 升 的 趋势 ,该 结论 
与 张志强 ”等 对 2000 一 2019 年 黄河 流域 植被 覆 羡 
度 时 空 变化 的 研究 结果 一 致 。 通 过 重心 迁移 模型 
发 现 ,各 级 林 草 植被 覆盖 度 重 心 都 呈现 出 向 西南 方 
向 迁移 的 态势 ,这 主要 是 因为 20 a 来 西南 地 区 的 地 
下 水 埋 深 与 年 均 温 的 变化 趋势 小 于 东北 地 区 ,降水 
量 的 变化 趋势 大 于 东北 地 区 (图 8)。 低 覆盖 类 林 草 
重心 在 整个 研究 时 段 均 位 于 永 济 灌 域 内 ,高 覆盖 类 
林 草 的 重心 主要 分 布 于 义 长 灌 域 的 西南 侧 , 各 级 林 
草 植 被 覆盖 重心 表现 为 从 东北 向 西南 方向 递减 。 
河套 灌区 2000—2019 年 林 草 生长 状况 整体 得 
到 改善 ,而 本 文 进一步 指出 了 在 这 20 a 间 河套 灌区 
林 草 地 的 数量 有 所 减少 ,这 依旧 是 恢复 河套 灌区 生 
态 环 境 的 较 大 挑战 ,同时 该 结果 也 与 王 娟 等 "得 
出 的 区 域内 植被 总 量 与 单位 植被 覆盖 度 并 没有 太 
大 的 相关 关系 的 结论 一 致 。 这 一 状况 一 方面 与 区 
域 自然 环境 关联 较 密 ,一 方面 也 与 区 域 社会 发 展 及 
工程 建设 相关 。 为 充分 利用 灌区 得 天 独 厚 的 条 件 ， 


人 


æ 0.02~0.04 mm 0.08~0.11 m 0.18~0.22 
me 0.04~0.06 æ 0.11~0.14 m 0.22~0.28 


年 降水 量变 化 趋势 
= -1.02~-0.66 ma 0~0.07 
 -0.66~-0.32 mm 0.07~0.16 m 0.53~0.78 


æ 0.31~0.53 


mm —0.32~0 m 0.16~0.31 mm 0.78~1.16 


0 25km 
L 


地 下 水 埋 深 变化 趋势 

= -0.14~0.01 m= 0.11~0.17 = 0.39~0.62 
= 0.01~0.06 mm 0.17~0.25 m 0.62~0.97 
mm 0.06~0.11 Æ 0.25~0.39 mm 0.97~1.46 


图 8 2000 一 2019 年 自然 因子 变化 趋势 
Fig. 8 Trends of natural factors from 2000 to 2019 
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自 建国 开始 ,灌区 开始 大 兴修 建 水 利 工 程 , 从 最 初 
的 “二 黄河 ”、“ 三 成 公 水 利 枢 纽 ” 到 现在 的 13 条 干 
渠 、339 条 支 渠 、 上 万 条 毛 汇 供水 至 灌区 ,使 灌区 水 
资源 得 到 了 更 均匀 的 配置 ,进而 也 有 利于 各 类 植 
被 的 生长 ,这 也 是 林 草 植被 生长 状况 逐渐 向 好 的 重 
要 原因 。 然 而 ,在 兴修 水 利 工 程 同 时 也 占用 了 大 量 
林 草 地 ,并 随 着 水 利 设施 的 完善 ,大 量 林 草 地 被 开 
F HAF , 20 a 间 新 增 耕 地 占用 林 草 地 718 km ,这 
也 是 导致 灌区 林 草 地 数量 逐年 减少 的 主要 原因 。 

本 研究 通过 地 理 探测 器 并 且 除 了 选取 部 分 自 
然 因素 ,还 选取 了 高 程 和 部 分 人 类 活动 数据 ,以 便 
更 全 面 性 、 针 对 性 的 分 析 影 响 林 草 FVC 空间 分 布 及 
变化 的 驱动 因素 。 研 究 发 现 , 距 水 域 距离 地 下 水 
埋 深 、 气 温 、 降 水 对 林 草 植被 覆盖 度 空 间 分 异 的 解 
释 力 要 相对 较 强 ,自然 条 件 是 影响 河套 灌区 植被 覆 
盖 度 的 主要 驱动 因子 ,而 水 热 条 件 是 影响 河套 灌区 
林 草 植被 生长 状况 的 主要 原因 ,这 也 与 李 成 林 等 3 
对 锡林郭勒 盟 干 旱 半 干旱 区 的 研究 结论 相似 。 并 
且 20a 来 , 随 着 气候 .地 下 水 等 因素 的 影响 ,西南 地 
区 的 较 高 植被 覆盖 度 和 高 植被 覆盖 度 的 数量 增加 ， 
而 东北 地 区 却 有 所 减少 ,通过 图 8 可 以 看 出 ,20 a 间 
西南 地 区 的 年 均 气 温 变 化 趋势 和 地 下 水 埋 深 变化 
趋势 总 体 小 于 东北 地 区 ,西南 地 区 的 年 降水 量变 化 
趋势 要 大 于 东北 地 区 ,说 明 土 壤 水 分 与 区 域 植 被 覆 
盖 度 呈正 相关 关系 ,而 气温 则 是 呈 负 相关 关系 ,这 
也 与 陈 效 还 等 3 对 内 蒙古 植被 覆盖 度 时 空 变化 的 
研究 结论 相似 。 人 类 活动 主要 是 距 居民 点 距离 . 距 
耕地 距离 和 距 道 路 距离 对 林 草 植被 覆盖 的 影响 较 
大 ,反映 出 人 类 活动 对 林 草 植被 的 影响 主要 体现 在 
开 虹 种 植 和 占用 建设 等 方面 。 

由 于 河套 灌区 主要 以 黄河 水 漫灌 为 主 ,因此 各 
灌 域 的 引 黄 水 量 与 土壤 盐 碱 化 程度 以 及 各 灌 域 的 
生态 治理 措施 都 有 着 较 大 的 差异 ,这 些 对 林 草 植被 
覆盖 度 有 着 重要 的 影响 。 但 在 数据 选取 方面 ,由 于 
数据 收集 的 限制 和 区 域 尺度 的 制约 ,本 文 未 能 加 入 
以 上 内 容 数据 ,以 上 影响 因素 依然 是 之 后 研究 需要 
重点 考虑 和 分 析 的 内 容 。 


4 结论 


=A 


通过 对 河套 灌区 林 草 面积 及 植被 覆盖 度 时 空 
分 布 特征 分 析 人 研究 ,得 出 如 下 结论 : 
(1) 2000 一 2019 年 河套 灌区 林 草 地 面积 总 体 呈 


减少 态势 ,减少 数量 为 966.15 km, 林 草 植被 覆盖 度 
等 级 整体 呈 上 升 的 趋势 ,多 年 平均 林 草 植被 覆盖 度 
为 28.3%。 空 间 上 ,河套 灌区 林 草 植被 覆盖 度 呈 现 
出 东北 向 西南 递减 的 空间 特征 ,变化 最 大 的 是 乌 兰 
布 和 灌 域 , 近 20 a 改 善 的 比重 较 大 。 

(2) 近 20 a 来 ,各 级 林 草 植被 覆盖 度 重 心 主要 
位 于 永 济 灌 域 东北 侧 和 义 长 灌 域 西南 侧 , 各 级 林 草 
植被 覆盖 重心 表现 为 从 东北 向 西南 方向 递减 ,并 均 
呈现 出 向 西南 方向 迁移 的 态势 ,各 时 段 迁 移 距 离 介 
于 0.71~15.46 km. 

(3) 距 水 域 距离 .地 下 水 埋 深 .气温 .降水 是 影 
响 区 域内 林 草 植被 覆盖 度 的 主导 环境 因子 ,对 林 草 
植被 覆盖 度 的 解释 力 均 大 于 0.3, 高 程 对 林 草 植被 覆 
盖 度 的 解释 力 最 弱 为 0.026; 各 驱动 因子 间 的 交互 作 
用 以 双 因 子 增强 为 主 , 其 中 距 水 域 距离 .地 下 水 埋 
深 、 年 平均 降水 与 其 他 因子 的 交互 作用 要 普遍 较 
高 ,解释 力 均 超过 0.35 ,说 明了 区 域 含水 量 的 大 小 是 
影响 河套 灌区 的 林 草 地 生长 状况 的 主要 影响 因子 ， 
通过 风险 探测 各 驱动 因子 对 河套 灌区 林 草 植被 生 
长 的 影响 均 有 其 适宜 的 范围 。 
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Spatiotemporal change and analysis of factors driving forest-grass vegetation 
coverage in Hetao Irrigation District based on geographical detector 


LI Xinlei', LI Ruiping’’, WANG Xiuging’, WANG Sinan', WANG Chengkun' 

(1. College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, Hohhot 010018, 
Inner Mongolia, China; 2. Inner Mongolia Key Laboratory of Big Data Research and Application in Agriculture and 
Animal Husbandry, Hohhot 010018, Inner Mongolia, China; 3. Surveying, Mapping and Geographic Information 
Center of Inner Mongolia, Hohhot 010050, Inner Mongolia, China) 


Abstract: Based on the data on forest and grass vegetation coverage from 2000 to 2019, the spatiotemporal 
change of forest and grass vegetation in Hetao Irrigation District over 20 years was analyzed, and the influence of 
various factors driving forest and grass vegetation coverage was quantitatively analyzed using a geographical 
detector. The results showed the following: (1) The area of forest and grassland in Hetao Irrigation District 
showed a decreasing trend in the studied 20 years, with a decrease of 966.15 km’. The vegetation coverage level 
of forest and grassland showed an overall upward trend. The average annual vegetation coverage of forest and 
grassland was 28.3%. On the spatial scale, the vegetation coverage of forest and grassland in Hetao Irrigation 
District showed a spatial variation characteristic of decreasing from northeast to southwest. (2) In the studied 20 
years, the centers of gravity of low forest and grass vegetation coverage and medium forest and grass vegetation 
coverage have mainly been located in the northeast of Yongji irrigation area and the southwest of Yichang 
irrigation area, respectively. The centers of gravity of forest and grass vegetation coverage at all levels show a 
decreasing trend from northeast to southwest, and all show a trend of migration to the southwest. The migration 
distance in each period was between 0.71 and 15.46 km. (3) Distance from water, groundwater depth, 
temperature, and precipitation are the dominant environmental factors affecting the forest and grass vegetation 
coverage in the region, with explanatory power of 0.427, 0.439, 0.318, and 0.368, respectively. The interaction 
between distance from water, groundwater depth, annual average precipitation, and other factors is generally 
higher, and the regional water content is the main driving force affecting the growth of forest and grassland in 
Hetao Irrigation District. 

Keywords: Hetao Irrigation District; vegetation coverage; spatiotemporal change; driving factor; geographical 
detector 


